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ABSTRACT 



We evaluate two dominant nuclear reaction rates and their uncertainties that affect 
production in explosive nucleosynthesis. Experimentally we develop thick-target yields for the 
'^°Ca(Q!,7)^'*Ti reaction at Ea = 4.13,4.54, and 5.36 MeV using 7-ray spectroscopy. At the 
highest beam energy, we also performed an activation measurement which agrees with the thick 
target result. From the measured yields a stellar reaction rate was developed that is smaller 
than current statistical-model calculations and recent experimental results, which would suggest 
lower ''''Ti production in scenarios for the a— rich freeze out. Special attention has been paid to 
assessing realistic uncertainties of stellar reaction rates produced from a combination of exper- 
imental and theoretical cross sections. With such methods, we also develop a re-evaluation of 
the '*^Ti(a,p)'*''V reaction rate. Using these two rates we carry out a sensitivity survey of ^''Ti 
synthesis in eight expansions representing peak temperature and density conditions drawn from 
a suite of recent supernova explosion models. Our results suggest that the current uncertainty in 
these two reaction rates could lead to as large an uncertainty in 44^j synthesis as that produced 
by different treatments of stellar physics. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — supernovae: general 



1. INTRODUCTION 



The dynamic synergy between observation, theory, and experiment developed over many years around 
the field of 7-ray astronomy has as its ultimate goal observations of specifc radionuclides informing our 
understanding of stellar explosions and the theoretical models that predict nucleosynthesis. Of the radioactive 
species observed so far, the most long lived, ^^Al and ^°Fe (with half-lives T1/2 = 7.17±0.24 x 10^ and 2.62 zfc 



0.04 X 10^ yr, respecti vely), are in reasonably good agreement with theoretical predictions ( Timmes et al 



1996; 



Diehl et al.l[2006h . Observations of those in the iron group, 56,57j^j (^^^^ ^ 6.075±0.01 d and 35.6±0.06 
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hr, respectively) and their decay products ^^'^"^Co (ti/2 = 77.233 ± 0.027 d and 271.74± 0.06d, respectively), 
are used in many ways to constrain our current models of the core collapse explosion mechanism. The 
radionuclide ^'^Ti = 58.9±0.3 yr), made in the same explosive environ ment but in much lower amounts 
compared to the very abundant nickle isotopes ( Wooslev et al. 1973 . 1995[) . is hoped to one day serve as an 
even more sensitive diagnostic and a valuable probe of the conditions extant in some of the deepest layers 
to be ejected. 

Observationally most of the attention has focused on the detection of the 68, 78, and 1157 keV 7-rays 
from the ''^Ti -> ^"'Sc — ^*Ca decay chain (see Figure 1). The 1157 keV 7-ray has been observed directly 
from a point source in Cassiopeia A (Cas A; lyudi n et al. 1994). T his was later confi rmed by observation 
of the low-energy ^''Sc 7-rays using the BeppoSax (jVink et al.ll200l[ ) and INTEGRAL (jRenaud et al.ll2006[ ) 
observatories. Using values for the distance, age, a nd 7-flux of Cas A, the amount of ^''Ti ejected was 
found to be 1.61°'^ xlO~^ Mq ( Renaud et al.l 120061 ) . in agreement with earlier observations using CGRO 
(jTimmes et al.lll996l ). Although the presence of ''"'Ti is currently below detection limits in SN1987A in the 
nearby Large Magellanic Cloud, its light curve is theorized to now be powered by the decay of ^''Ti. The 
yield of ''^Ti in SN1987A has been estimated from its light curve to be 1 — 2 x 10~^Mq, a factor of 3 greater 



than predicted by models ( Diehl et al. 20061) . A third '^'^Ti source ( of lower significance compared to the one 
in Cas A) has been observed in the Vela region (jivudin et al.lll998l ) but the existence of a co- located young 
supernova (SN) remnant has not been confirmed. While the mass of 44^j observed in SN remnants appears 
to be underproduced by past models, the number of observed sources of ■*'*Ti in all-sky surveys appears to 
be less than expected from estimates of the Galactic S N rate and the known '*'*Ti half life, leading some to 
question whether '^''Ti-producing SNe are exceptional (jThe et al.ll2006l ). 



Theoretically, '***Ti production is traditionally ascribed to regions experiencing a strong "a— rich freeze 
out" , where material initially in nuclear statistical equilibrium (NSE) at relatively low density is cooled so 
rapidly that free a— particles do not have time to reassemble, through inefficient 3— body reactions that span 
the mass gaps at A=5 and A=8, back into the iron group. A likely scenario is material in or near the silicon 
shell as it experiences shock wave passage during a core collapse SN event. Assuming the material cools 
adiabatically over a hydrodynamic timescale (thd = 446x/^//o7), where x is a scaling parameter (here unity) 
and Pi is the initial (peak) density in g cm^'^, initial conditions that would result in a final a-particle mass 
fraction of ^ 1% would require temperatures high enough to ensu re NSE {Tgi ^ 5) and an initial density 

< min(4.5 x lO^Tg^, 2.5 x lO^Tg^x^ 
greater than unity (Section 3.2). 



) (IWooslev et al.lll973l ). This translates to a radiation entropy 



In models of massive stars the a— rich freeze out dominates the solar production of several species, includ- 
ing ''^Ca (made as Ti), "^^Sc, ^^Fe (made as Ni), and ^^'^^Ni, while still others seem to r equire a sizable com- 



ponent to account for their solar abundanc es, including ^Cr, ^^Co, °^Ni, and °^Zn (jWooslev et al 



199£ 



Thielemann. Nomoto. fc Hashimoto 1 119961 ) . Although ^^Ni is the dominant species produced in both NSE 



and the a— rich freez out over a wide range (0 < ry < 0.01) of neutron excess (jHartmann. Wooslev. fc El Eid 



19851), the solar abundance of ^^Fe is don iinated by production in explosive s ilicon burning in massive stars 
and by a large contribution from SNe la (jTimmes. Wooslev. fc Weaveiill995[ ). 



In order to address one aspect of the model uncertainties associated with the theoretical predictions of 
"^^Ti in SNe, we focus on exploring the nuclear data uncertainties of two key reaction rates. Experimentally 
we address the ^°Ca(a,7)'*^Ti cross section where the existing experimental results are inconsistent and 
theoretical estimates are co mplicated by the sup pression of El T = — >■ T = gamma transitions in self- 
conjugate {N — Z) nuclei ([Rauscher et al.l 120001) . This cross section has been measured in the past using 
several techniques. Cooperman et al. |l977l ) "used in-beam 7-ray spectroscopy to measure the capture of 
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a-particles on a metallic calcium target in the centcr-of-mass energy range Eqm — 2.5 — 3.65 MeV. In that 
work the excitation function for the 1083 keV firs t excit ed state transition in ''^Ti was determined and the 
resonance strengths were developed. Nassar at al. ( 20061 ) determined the integral cross section in the range 
EcM — 2.2 — 4.17 MeV by bombarding a He gas target with a ^°Ca beam and collecting the recoiling ''^Ti 
in a catcher foil. Accelerator mass spectroscopy was then used to determine the ratio of '*'*Ti/Ti from the 
known content of Ti in the catcher foil. Recently, a slightly broader energy range of E cm = 2.11 — 4.19 MeV 
was explored using the DRAGON recoil mass spectrometer (IVockenhuber et al.ll2007l ) in inverse kinematics. 
They developed individual resonances whose sum was used to determine a reaction rate for 1 < Tg < 5.5. Yet 
despite these often heroic expenditures of time, toil, and treasure, in the temperature range of astrophysical 
interest there still exists a factor of 3 or more difference between the experimentally determined reaction 
rates. 

In this work, we develop the cross section for '*"Ca(Q!,7)'''*Ti by two separate methods as a check on 
systematic uncertainties. First we used in-beam 7-ray spectroscopy to measure a thick target yield of the 
^''Ca(a,7)*'*Ti reaction. We then determined the number of '^'^Ti nuclei produced by counting low-energy 
7-rays from the decay of ^"^Ti in an irradiated target. Special attention was devoted to checking the internal 
consistency of the measurements and to establishing realistic uncertainties in developing the stellar reaction 
rate from a combination of experimental and theoretical cross sections. We have made a similar evaluation 
of the stellar react ion rate for the domin ant destruction reaction, '*'^Ti(a,p)^^V, based on the original ex- 
perim ental work of ISonzogni et al. and the theoretical cross section work of iRauscher fc ThielemannI 
(|200lh . 



Our results are presented in four parts. In Section 2, we describe our experimental efforts. Section 3 
then discusses the development of stellar reaction rates, and an expose of the past and present experimental 
and theory efforts for the two rates in question. In Section 4, we present nucleosynthesis results for the 
production of ^'^Ti and 56.57,58]\jj reported in previous surveys of massive star evolution. We then carry out 
a sensitivity survey of ''^Ti to variations in the principle production and destruction rates using more recent 
SN models. In Section 5 we provide a discussion and conclusions. 



2. Experimental Methods 

The '*°Ca(a,7)'*''Ti reaction was measured at the Lawrence Livermore National Laboratory (LLNL) 
Center for Accelerator Mass Spectroscopy (CAMS) using a 10 MV FN Tandem Van de Graaff accelerator. 
To calibrate the beam energy a silicon detector was placed in a low current beam and the measured energy 
was then compared with spectroscopy grade a-sources. The a energies used were the 5.30 MeV from ^^"^Po, 
the 6.12 MeV from 252cf, the 5.49 MeV from 24iAm, and the 4.69 and 4.61 MeV a's from the decay of 
^^°Th. From the calibration, the uncertainty in the a-beam energy is ± 5 keV. The thick target yield was 
measured at Ea = 4.13, 4.54, and 5.36 MeV beam energies. 

For each beam energy a target was manufactured by pressing "^'CaO powdei0 into a copper holder. The 
powder had a purity of 99.95% (metals basis) but contained ppm concentrations of C and F. To completely 
stop the beam, each target had a minimum thickness of at least 1.1 mm. The target was mounted on a 
copper block within the vacuum chamber and tilted at 30° with respect to the beam. The vacuum chamber 
contained two windows so that the target could be visually monitored. A schematic of the target chamber 



1 Obtained from Alfa Aesar, MA, USA 
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is given in Figure [2] 

The target chamber was electrically isolated from the rest of the beam line. The target was electrically 
connected to the target chamber allowing the beam current to be directly measured from the chamber. 
The current integration from the target chamber was checked to a precision of better than 1% using a 
NIST-traceable precision DC current source. An opposing-pair magnet was attached upstream of the target 
chamber to suppress the escape of secondary electrons generated by the beam on the target. A summary of 
the irradiation runs is given in Table [TJ 

Two 80% HPGe detectors were used to measure the prompt 7-ray yield during irradiation. One detector 
was 11.4 cm away from the target at 30° with respect to the beam. The location of the second detector was 
at 15.8 cm and an angle of 99°. The HPGe detector thresholds were set at 125 keV. The average deadtime 
during a run was 10 %. 

The 7-ray energy spectra was accumulated in 8192 channels using two Ortec AD413a ADCs. Efficiency 
and energy calibrations of the HPGe detectors were made using ^°Co, ^^Na, ^^^Cs, ^''Mn, and ^'^•^Ba NIST- 
traceable sources with activities known to a la uncertainty of 1%. Efficiencies of 0.085% and 0.162% were 
found for the detectors at 30° and 99°, respectively. 



2.1. Analysis of Prompt 7-Ray Data 

The thick target yield for the *°Ca(a,7)*''^Ti reaction was deduced from the yield of the 2+— )-0"'" 1083 
keV transition which collects most of the strength from ''^'^Ti. The angular distribution of the 1083 keV 
7-rays with respect to the angle between the 7-rays and beam direction is given by W{6) = J^,, ai Pi (cos 



with I = 0, 2, 4. The cross section is proportional to oq so only this term needs to be determined (jPver et al 



19811 ). By placing the detectors at angles for which P4(cos0) is zero, 30.6° and 109.9°, uq can be determined 
from measurements at only two angles. The experimental thick target yield of the 1083 keV 7-ray at angle 
is given by 

Y^os3{e) = ^ , , (1) 

IVpLtAi, £1083 

where Nc is the number of counts in the 1083 keV photopeak, Np is the number of a particles impinging 
on the target, Lt is the detector live-time fraction. At is the natural abundance of ^''Ca, and eioss is the 
efficiency of the HPGe detector at 1083 keV. The region near the 1083 keV 7-ray is shown in Figure[3]for Ea 
= 5.36 MeV. The 1083 keV photopeak lies on the tail of the 1039 keV ^"Ge doppler shifted 7-ray excited by 
fast neutrons primarily from the ^^F(a, n) reaction. The bac kground and 1083 keV peak were fit to an error 



function convoluted with a Gaussian, as in lGete et al.l ()l997l ). The total experimental yield for the 1083 keV 



7-ray was determined by finding ag from the angular distribution of Y{0). The experimental yield of the 
1083 keV 7-ray determined from these fits are given in Table [H In order to convert the yield of the 1083 keV 
7-ray into the yield of ''^Ti, one must take into account those transitions which bypass the 2+ 0+ 1083 



Table 1. Irradiation Conditions for the "'^'CaO Targets. 



Irradiation Run 


Bo (McV) 


Irradiation Time (hr) 


Total Charge (^iC ) 


1 


4.13 


6 


45808.6 


2 


4.54 


7 


46115.8 


3 


5.36 


10 


25763.2 
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keV transition in ''^Ti. The Monte Carlo program DICEBOX (jBecvaiill998l ) was used to simulate the 7-ray 
cascades from the decay of the compound nucleus 44^j jj^ order to estimate the number of transitions which 
bypass the 2+ — ?> 0+ transition. These simulations suggest (20 ±3)% of transitions bypass the 2+ — > 0+ 
transition with the error on the simulations taking into account the uncertainty in the nuclear level density, 
the photon strength function, and the capture state. In Figure |4l the thick target yield corrected for the 
missed strength to the 1083 keV level is compared to the NON-SMOKER cross section. 



The experimental yield can be related to a theoretical cross section a{E) by the equation 

-dE/dx 



(2) 



where ct{E) is the energy-dependent cross section, and dE/d x is the stopping p ower for "^*CaO. The dE/dx 
values for "'^'CaO were calculated using the program SRIM (jj.F. Ziegleiil2004l ). Table [2] gives a comparison 
of the e xperimental and calculated yie lds using the NON-SMOKER Hauser-Feshbach '*°Ca(Q!,7)'*'*Ti cross 
section faauscher fc Thielemann 2001 ). The experimental yields are a factor of 1.7-2.3 times smaller than 
the yield calculated from the theoretical cross sections. 

The uncertainties from detector efficiencies, deadtime, angular distribution corrections, coincidence 
summing, and beam current integration are tabulated in Table [3l The integrated beam current was cross 
checked by using the ^'^0{a,a'j) and ^"'Ca(Q;, q;'7) reactions. The thick target yie ld for the ^''0(a,Q;'7) 
reaction that excites the first excited state at 870.8 keV was measured in 'HsuI |l982[) at E^ = 5.486 MeV. 



Using Equation (1), Np from ^'^Of a, a'"| ') was determined to be (7.2±0.86)x 10"^^ after correcting for the 



HsuI 



difference in beam energy between iHsul (|1982I ) and our work. The thick target yield for the Coulomb 
excitation of the 1157 keV level in ^''Ca was calculated using the program GOSIA ijCzosnyka et al.l 119831 ) 
and used to infer Np of (8.2±1.6)x 10^^. This compares well with our beam current integrator which gave a 
value of 7.2x10^® for Np. The efficiency of the detector at 30° was sensitive to the decay location due to the 
attenuation of 7-rays through the copper target holder. By varying the position of 7-ray sources to match 
the approximate 1 cm^ beam spot size the uncertainty in the efficiency was determined. The uncertainty due 
to coincidence summing of 7-rays in a single detector was estimated from the geometric solid angle spanned 
by the detectors. 

The angular distribution is the largest source of uncertainty. Our detectors did not sit at exactly the 
zeros of the P4 term in the angular dist ribution which introd uces some uncertainty when aq is determined 
with only two detectors. Furthermore, ISimpson et al.l (jl97ll ) measured the angular distribution of 7-rays 
following resonance a capture at Ea = 4.22, 4.26, and 4.52 MeV and found a strong contribution for the 
P4 term. The uncertainty on ao was estimated by the following procedure. A value was chosen for the 04 
coefficient between ±0.5 and then a fit of the angular 7-ray yield was made. The resulting variation in oq 



Table 2. '"^Ca(Q;,7)'''*Ti Measured and Theoretical Thick Target Yields per a Particle on Target. Yioss is 
the measured yield for the 1083 keV 7-ray, YnTi is the total yield for the productio n of "^^Ti, and Ithnorv the 
thic k targe t yield calculated from Equation [2] assuming the theory cross section of lRauscher fc Thielemann 
(|200lr ). Fofflino is the yield from our activation measurement of the target irradiated at 5.36 MeV. 



(McV) 


Yl083 (10-") 


Ki4Ti (10-") 


^theory (10-") 


yottlinc (10-") 


4.13 


2.11 ± 0.40 


2.53 ± 0.50 


5.96 




4.54 


5.72 ± 1.1 


6.86 ± 1.4 


16.4 




5.36 


29.2 ± 5.7 


35. 0± 7.1 


61.0 


35.7 ± 2.5 
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was between 17% and 18% for different values of a^. This is a conservative estimate because the many 7-ray 
cascade paths and alpha energies involved are likely to wash out any nuclear alignment and would cause 04 
to be nearly zero. 



2.2. Low Background Counting 

The offline counting of the irradiated target took place at the Low Background Counting facility at 
LLNL. Only the target irradiated at Ea = 5.36 MeV was counted because the activity of the other targets 
was estimated to be too low. A HPGe low-energy photon spectrometer (LEPS) detector was used to detect 
68 and 78 keV 7-rays from '*''Ti decay. The target was placed 2 mm away from the detector face and counted 
for two weeks. The activity of the ^'^Ti target was determined by comparing it to a 56.6 ± 1.6 nCi ^^Ti 
reference source. The reference source was counted with the LEPS detector in the same geometry as the 
^^Ti target. 

The ^''Ti reference source strength was determined by using an 80% HPGe detector to compare the 
1157 keV 7-ray to a ^^Na and ^°Co calibrated source. The count was made at a distance of 60 cm away 
from the detector in order to avoid summing of the 1157 keV 7-ray with the 511 keV 7-ray. 

The yield was found from 

ATi 

Y = 2-, (3) 

7Vpln2' ^ ' 



wher e A is the activity of the irradiated target and T1/2 — 58.9 ± 0.3 years is the Ti half-life (jAhmad et al 



20061 ) 



A simultaneous fit of the peaks in Figure [5] gives an experimental yield of (35.7 ± 2.5) x 10^^^ ''^Ti per 
a-particle. The uncertainty in the offline yield takes into account the uncertainty in the detector efflciency, 
half-life, integration of beam current, calibration source strength, and statistics. Since the activity of the 
'*'*Ti target was determined with a reference source having the same decay scheme the need to correct for 
summing of the 68 and 78 keV 7-rays was eliminated. This yield is 22% higher than the yield from the online 
counting and agrees with the estimate that (20±3)% of 7-ray cascades bypass the 1083 keV level. This also 
demonstrates that sputtering of the target was minimal during particle bombardment. 

Table 3. Compilation of systematic uncertainties for the online measurement (Icr). 



Source of uncertainty 


Uncertainty 


Beam integration 


1% 


Detector efficiency 


6%-7% 


Deadtime 


1% 


Angular distribution 


17%-18% 


Coincidence summing 


1% 
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3. Reaction Rates 

3.1. The True Gamow Window 

Thermonuclear reaction rates are obtained through inte gration of the energy-dependent cross section 
weighted by a Maxwell-Boltzmann (MB) particle distribution (|Fowler. Caughlan fc Zimmerman .ilQGTi) . For 
reaction /^(j, k)L'^ 

3.732 X IQi" 



A 



1/2 



rp 3/2 



a'^^E'; exp{-lim5E''/Tg)dE[ 



(4) 



where Na is Avagadro's number, Aj ^ 3.64 is the reduced mass of the target and incident projectile in 
atomic mass units, n is an index representing bound states in the target and product nuclei (/i = for the 
ground state, 1 for first excited, etc.), and v is the relative velocity of the a-particle (j) and the target (/^) 
in cm s~^. The cross section tr^^ is expressed in barns, i^j* is the center of mass energy of the a-particle and 
the target in MeV, Tg is the temperature in billions of degrees Kelvin, and the reaction rate has units of cm^ 



There are two important issues regarding the cross section a^f. used in the integration and its relation to 
the experimental cross section. First, when computing astrophysically relevant reaction rates, the cross sec- 
tion has to be the one for a thermally excited target nucleus in t he plasma. Only such a stellar cross sectio n 
allows application of detailed balance to derive the reverse rate ( Holmes et al. 1976 : Rauscher et al. 20091) . 
Laboratory cross sections are measured with the target nucleus in the ground state only and the correction, or 
stellar enhancement factor, has to be calculated from theory. For the ^°Ca(a, ^Y'^ Ti and ^"^Ti(a, p)^^V reac- 
tions the correction factors are unity over the entire temperature range of interest (| Rauscher fc Thielemann 



20001 ). Therefore the cross section cr^^ is identical to the cross section a appearing in Eq. [2l 



Second, the question arises as to how well the reaction rate is constrained by the experimental cross 
section data. In the rare case of a thick target yield measurement that completely samples the relevant stellar 
temperature range, a direct calculation of the thermonucl ear reaction rate can be obtained by introducing 
a simple change of variable and combining Eq.'s [2] and |4] (IRoughton et al.lll976l ). More often limitations of 
the experimental apparatus and/or beam time allotment restricts the measurement of cross section data to 
a fairly narrow energy range. To produce a reaction rate spanning a realistic temperature range for stellar 
synthesis, one often has to supplement the experimental results with theory cross section values to perform 
the integration in Eq. 31 The integrand exhibits a maximum contributing most to the reaction rate due 
to the folding of the energy-dependences of the cross section and the MB distribution. Conventionally, the 
location and width of this Gamow peak are estimated by multiplying a Coulomb barrier penetration factor 
with the high-energy tail o f the M B distribution or by considering a Gaussian approximation to it (see, e.g., 
ClavtonI |l984l ): lRauscherl I2OI0I )). 



However, this is applic able for capture reactions only as long as the 7-width is larger than the charged 
particle width (|lliadisll2007l ). Since the a width is changing rapidly with energy due to the Coulomb barrier, 
this prerequisite may only be fulfilled at very low projectile energies which effectively shifts the Gamow 
window to lower energies than expected from the standard approximation formula. Figure [6] confirms this 
by showing the "true" Gamow window at Tg = 2 and Tg = 5 in comparison to the peaks obtained from 
the standard formula and its approximation for the ^'^Ca(a, 7)^^Ti and ^^Ti(a,p)^^V reactions. Finally, 
Fig. [7] compares the energy range spanned by our experimental data with the true Gamow window for 
stellar temperatures 1 < Tg < 5. The width of the window is defined such that filling the "window" with 
experimental data would determine the reaction rate to 10% accuracy. The minimum temperature at which 
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an experiment can fully inform the reaction rate is the overlap of the full width of the window and the range 
of experimental data (shaded region). For the astrophysical processes considered here, the most relevant 
temperature range is between 2.0 < Tg < 4.0 (Section 3.4). A less severe situation exists for the '*'*Ti(Q;,p)^^V 
(Figure [5]). From these figures it is evident that the limited range of experimental data collected in these two 
experiments do not contribute significantly to the rates in the important temperature range and that the 
supplemental theory cross sections at lower energies will dominate both of the rates and their uncertainties. 
Nevertheless, the data can be used to test the NON-SMOKER predictions at higher energies. 



3.2. Reaction Rate for 40Ca(a, 7)^4X1 



To obtain a "semi-experimental" astrophysical rate for '^'^Ca{a, 7)"'^Ti, we normalized the NON-SMOKER 
cross section by a factor of 1.71 to agree with our measured thick target yield at 5.36 MeV (Figure |4|). The 
data at this energy were chosen because the experimental result was confirmed by two independent tech- 
niques. Although the errors on the experimental thick target yield are small, during the calcuation of the 
reaction rate we have to assign an uncertainty which is similar to the uncertainty on the theory cross section. 
Therefore we will perform our astrophysical studies with an assigned factor of 2 uncertainty on this rate. 
This also accounts for a potentially different energy-dependence of the theoretical cross section, and fits well 
with the general factor of 2 — 3 accuracy exp ected for global predictions of low-energy a-capture reactions 
on self-conjugate nuclei (jRauscher et al.ll2000l ). 



The resulting reaction rates were then fit assuming the REACLIB parameterization (jRauscher fc Thielemann 



20011) 



Na < cfv >= cxp(ao + aiTi 



9 



rj, -1/3 



+ asTg^/^ + a^Tg + a^Tg^^^ + ag 111(^9)) 



(5) 

where Tg is the temperature in 10^ K. The fit is defined in the range 0.01 < Tg < 10.0. The reverse rate for 
^'*Ti(7, Q!)'*''Ca is calculated through detailed balance, the Q value of the forward reaction is 5.1271 MeV. 

Figure [9] shows the ^'^Ca(Q!, 7)'*'*Ti reaction rates considered in this study. The values are generated 
from fits to reaction rates developed from numerous experimental and theory efforts carrie d out over the 



last 2 years. Some have been used i n widely cited n ucleos ynthesi s surveys. In particular. IWooslev et al 



(|l995l ) utihz ed the theory rates froni IWoosley et al.l (119781). while Limongi fc Chieffi I (j2006r ) used the fit 



published in Rauscher fc Thielemann ( 2000l ). Rauscher et al. 



Tg = 3 to the " empirical" rate developed in iRauscher et al 




used this same fit but normalized at 
. Only "experimentally" determined 
ra tes are shown in the main pan el, along with the tabulated Hauser-Feshbach theory based rate (crosses) 
of iRauscher fc ThielemannI (j200l[ ). hereafter referred to as NON-SMOKER. Also plotted in the inset is the 
ratio of each rate versus the NON-SMOKER rate, our recommended '*''Ca(Q;, 7)'*^Ti rate lies between the 
oldest and most recent experimental results. 

For ^''Ti synthesis the most relevant temperature range is 2.0 < Tg < 4.0. Over this range experimental 
'^°Ca( o;,7)^"^Ti rates vary by fac tors of 3-5, respectively. The two most recent experimental rates (jNassar et al 



2006 



Vockenhuber et al.l 120071 ) are larger than our recommended rate an d would suggest increas ed Ti 



synthesis over our result and the yet smaller original experimental result of ICooperman et al.l (|l977n . 



The range of variation in the theory rates is also large, tha t of Wooslev et al. ( 1978 ) being of the order of 
the largest experimental results, the smallest is clearly that of Rauscher fc Thielemann ( 2000l) . Most unfor- 
tunately, the rate generated from this published fit differs from the tabulated NON -SMOKER rate (crosses) 
by fac tors of 2 (high and low) on both ends of the important temperature range. IRauscher fc Thielemann 



(|2000l ) fit this rate to an accuracy of 0.92 over the temperature range 0.1 < Tg < 10. In general fits to 
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charged particle rates have larger deviations than those for neutron-induced reactions. Most often a quote of 
low accuracy pertains to deviation at the lowest temperature points of the fit. Base d on argumen ts related to 
the nu clear level density and its impact on the applicability of the statistical model, Rauscher fc Thielemann! 



(120001) suggest that the lower limit of applicability for this dominant production rate is Tg = 0.24, well below 
the important temperature range in this study. To be fair the statistical result itself could well be in error 
by a factor of 2 or more due to the uncertainty in the global a+nucleus optical potential, and the fit does 
agree at Tg = 2.8, which is close to the midway point of ''^Ti production in several of the scenarios studied 
in the next section. We have re- fit this rate (dotted line in the inset), and observe that over the important 
temperature range it now differs from the NON-SMOKER rate by no more that 12% (typically better than 
5%). We will include both fits in our sensitivity analysis. 

With an uncertainty of a factor of 2, our new "semi-experimental" rate encompasses most of the pre- 
vious experiments and calculations. This underlines the necessity for further measurements in the relevant 
astrophysical energy range (2-5 MeV) . 



3.3. Reaction Rate for ^'^Ti{a,p)'^'^V 



Figure fTOl s hows a similar comparison for t he '*''Ti(a,p)'*^V rates. As with the dominant production 
rate, the original Rauscher fc Thielemann ( 2000l ) fit to this destruction rate is in poor agreement with the 
statistical NON-SMOKER rate (crosses) . They list a fit uncertainty of 0. 12 and the lower limit of applicability 
as Tg = 0.035, again outside the important temperature range in this s tudy. We also note th at the NON- 



SMOKER rate exceeds the older (dot-long-dash line) statistical result of lWooslev et al.l (11978 ) for Tg > 2.4, 



which would suggest decreased Ti production over that seen in the survey of lWooslev et al 



(199 



i. 



The only published experimental reaction rate (jSonzogni et al.ll2000l ) is roughly a factor of 2 above the 
NON-SMOKER rate, which would suggest even further decreased '^^Ti synthes is compared to the previously 
available theory rates. To generate their reaction rate Sonzogni et al. ( 2000l ) interpolated their four cross 
section points between Ec.m. = 5.7 and 9. MeV, extended the energy range abo ve and below these values 
with scaled cross sections from SMOKER (jThielemann. Arnould. fc Truran 1119871 ). and then integrated over 
a MB distribution to obtain the reaction rate and produced a REACLIB fit to it. They claim agreement 
with the calculated rate of 43% at the highest temperatures and 30% at Tg — 2.5. Interestingly, their 
de rived experimental rate s hows a ve ry different temperatu re dependence compared to the reaction rates 
of Rauscher fc Thielemann (2001) and Wooslev et al. ( 19781 ). suggesting a very different energy dependence 
of the cross section. This result is very surprising as the ener gy dependence of the S MOKER cross section 
is very similar to that of the NON-SMOKER cross section. Sonzogni et al. ( 2000[) provide plots of their 
experimental values and calculated c ross s ection between 5 and 10 MeV, and the fit of their reaction rate 
vers us SMOKER an d Wooslev et al.l (jl978l ). but unfortunately no further details. Considering Figures [6] and 
|8] of [Sonzogni et al.l hoOfll ). it is clear that their reaction rate for 2.5 < Tg < 5 was determined completely 



by the scaled SMOKER cross section values. 



We provide a re-evaluation of the rate from the data of ISonzogni et al.l (|2000( ) which also includes an 
increased uncertainty estimate. Since the experimental data (within errors) lies within 20% of the NON- 
SMOKER cross section prediction, we generated a set of hybrid cross sections by supplementing the exper- 
imental ones with unrenormalized NON-SMOKER values above and below the experimental energy range. 
We assigned a factor of 2 uncertainty to the theory values above the experimental data and a factor of 3 
below. For the experimental data points we assigned the errors on energy and cross section as given in 
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SonzQgni etall |2000l) . The code EXP2RATE |Rauscheij|2003-2009|) allows one to calculate reaction rates 
from cross sections (or astrophysical S'-factors) that include both theoretical uncertainties and experimental 
error bars. The ratio of the resulting rate with its uncertainty range to the NON-SMOKER rate is shown 
in Fig. 111! Again, the experimental data contribute to the rate integral only at higher temperatures and 
the standard NON-SMOKER rate is still well within the uncertainty even at those temperatures. In our 
sensitivity study we adopt the geometric mean of the upper and lower rate limit as our "recommended" 
rate and consider a variation of a factor of 3 (up and down) as a reasonable uncertainty estimate based on 
inspection of the uncertainty limits in the relevant temperature range. 

Note that our semi- experimental rate for ^'^Ti(ct, p)'^^V is fit in the en doergic {—Q) direction while 
the original theory rates (jRauscher fc ThielemanrJl2001t IWooslev et al.lll978[ ) were fit in the exoergic {+Q) 
direction. Reaction rates for targets in their ground states alone (or any distribution other than the rmal equi 



libriu m) do not obey reciprocity because the forward and reverse reactions are not symmetrical (H olmes et al 



19761) . Thus it is very important to measure cross sections for astrophysical application in the direction that 
is least affected by excited state effects in the target, which is almost always in the exoergic direction. This 
particular reaction is an exception in that it has a small reaction Q value (—0.410 Me V) and excite d state 
effe cts that are minimized by Coulomb suppression of the stellar enhancement factor. See lKiss et. al. I (|2008. ) 
and Rauscher et al. ( 20091 ) for details. 



In Table JH we provide fits in the REACLIB format for all of the rates discussed above. For nasr06 
(|Nassar et al.ll2006[ l '.he reaction rate values we fit were themselves generated from a fit supplied in that paper, 
otherwise we have used (a,nd alw ays recom mend) tabulated yalues where available , as we have done for rathOl 
iRauscher fc ThielemannlboOli l. coop77 (jCooperman et al.l[l977l ). and wfhz78 (jWooslev et all 1 19781 ). The 
rates lab eled "hsrlO" are our ev aluated rates from our own experimen tal work and our re-evaluation of the 
rate from lSonzogni et al.l (|200Cll ). For voch07 (jVockenhuber et al.ll2007l ) we fit their central "Rate" tabulated 
between 1 < Tg < 5.5 (see their Table III). For each fit we define a measure of accuracy, C, between the rate 
Ti and the fit fi for the seven Tg temperature points 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, and 4.0 as 



c 



i=l 



(6) 



We also supply the reverse rate fits in the REACLIB format. Thes e must be multiplied by t he ap - 
propriate ratio of nuclear partition functions as specified and tabulated in IRauscher fc Thielemann! (|2000f ) . 
Again, the re-evaluated and re- fit rates for ^''Ti(Q!,p)^^V must be inserted into REACLIB in the forward 
rate direction (increasing mass or charge). 

Of course m any other nuclear reaction rates affect "^^Ti synthesis. A prioritized list was suggested by 
The et al. ( 19981 ). but only for one choice of peak temperature and density (Tgp = 5.5, pp = 1.0 x 10^ g 



cm~'^) and three values of electron fraction (Yg = 0.5,0.499, and 0.497). Recently the NuGRID team has 
presented a preliminary survey of ^^Ti synthesis over a wide range of peak temperature and density conditions 
for Ye = 0.5 in which they define several regions where a similar set of reaction rates are suggested to be 
more effective than others (jMagkotsios et al.1120081 ). 



For this study we used the reaction rate library developed bv iHoffman et al.l (|2002[ ). We only explore 



variations in ^'^Ti synthesis due to the dominant product ion and destructio n ra tes mentioned above, but here 



cite our sources for several other key reactions noted byEhl^ & and kagkotsios et al.l & 



For important rates that produce ''^Ti we adopted the 3a reaction rate of ICaughlan fc Fowled (|l998l ). 
Coupled with their ^^C(a, 7)^^0 rate uniformly multiplied by a factor of 1.7 (corresponding to an 5- factor at 
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300 k eV of 170 keV barns), these proved optimal for producing the solar abundance set (|Timnies. Wooslev. fc Weaver 
19951). Although we do not explore variations in these rates, we acknowledge their importance, especially 
the former, in setting the a-abundance on which both '^'^Ca and '^^Ti production depend. This choice is not 
critical to showing the sensitivity of "^^Ti to the two rates we concentrate on. As previously mentioned we 
developed a new semi-experimental rate for ■**^Ca(Q;,7)'*^Ti with a factor of 2 uncertainty and also consider six 
alternate reaction rates (Figure[9l). We also adopted rate s for the weakly competing side chain ^°Ca(a,p)^'^Sc 
and '^^Sc{p,j)^'^Ti from iRauscher fc ThielemannI |2000l ). 



The dominant "^^Ti destruc tion rate in this study was *'*Ti(a,p)'*^V. We explore three choices, our re- 
evaluated Sonzogni et al. ( 2000l) rate with uncertainty factors 3 and 1/3 as upper and lower limits . Othe r 
destruction reactions of note were '^'^ Tifp, 7)^^V and "^^VfiP, 7) ^ ^Cr, both taken from IPisker et all ( 2001 ). 
Finally *'^Ti{a,-f)'^^Cj: was taken from IRauscher fc ThielemannI |2000l ). 



In the following section we adopt our recommended "semi-experimental" rates as the default production 
and destruction rates when we report specific values in our nucleosynthesis survey. 



4. Nucleosynthesis Studies 

4.1. Initial Conditions 

Previous efforts to explore the sensitivity of '^''Ti synthesis to the ^°Ca(Q!,7)^'*Ti reaction rate have 
often been confined to one-zone parameter ized network simulations for a limited number of initial conditions 



(|The et al.l 119981 : IVockenhuber et al.l 120071). A much larger survey promises to explore Ti synthesis over 



a wider range of initial conditions (jMagkotsios et al.l 120081 ). From this survey we choose eight points that 
reflect peak conditions experienced in various SN explosion models that span the regimes from incomplete 
silicon burning to the a-rich freeze out (see Table [5l). Ideally one would prefer to explore rate sensitivity in 
full star models of stellar evolution ( Nassar et al. 20061 ) . but for species that are made under a limited range 
of stellar conditions that do not involve details affected strongly by the stellar physics, such as burning over 
long time-scales in convective regions, one zone calculations are usually adequate as long as a reasonable 
range of conditions are explored. We take this approach. 

Our calculations are straight-forward. Table [5] lists our choices for peak temperature and density. 
From these we calculate the hydrodynamic time-scale, thd — 446x/-y/p^ and the radiation entropy, S^ad = 
3.33 X T^pl P5p, where p^p = 10~^ 'x Pp g cm~^ and S is in units of Boltzmann constant per baryon (hereafter 
we refer to entropy without its units) . For the hydrodynamic timescale x is a scaling parameter (assumed to 
be unity unless explicitly stated). The material then expands adiabatically (p oc T'^) with p{t) — ppexp(— i/r) 
on a hydrodynamic time-scale until the temperature declines to Tg 0.25, a point where all charged particle 
reactions affecting ^'^Ti synthesis have frozen out. The last column (t^) gives the approximate time for each 
simulation to terminate based on the scale factor for the hydrodynamic time-scale. We will explore two 
values, X = 1 and 5. The initial compositions consist of neutron, proton, and a— particle mass fractions that 
provide for a specific neutron excess rj = ~ Zi){Xi/Ai) where Ni, Zi, Ai and Xi are the neutron, 

proton, atomic mass number, and mass fraction of the isotope i with "^^Xi — 1. The neutron excess is 
related to the electron mole number via Yg = 1 — 2ri. For each initial condition we will survey a range of 
neutron excess (—0.01 < rj < 0.03) as well as variations in two specific reaction rates. 

Our peak conditions were chosen t o reflect the sensitivity of the nucleosynthesis to variations in entropy 
and expansion timescale. According to Magkotsios et al. ( 20081 ). points CasA 1 — 3 lie along the track of a 
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Table 4. Reaction Rate Fit Parameters for ^^^C&ia, 'Y)^^Ti and ^^Ti(a,p)^^V 
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Table 5. Nucleosynthesis Survey - Peak Initial Conditions 
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model for Cassiopia A and have very similar entropy values near 30, they differ from one another by a factor 
of 1.5 (33%). Points 2DExpl 4 — 6 are from a rotating two-dimensional explosion model (point 5 closely 
corresponds to a point along the track of a Gamma-ray Bm^st model), while points 2DMHD-7 and -8 lie along 
the track of a two-dimensional rotating MHD star (point 7 samples the "Chasm" region where '*'*Ti synthesis 
is suppressed). We chose points along these tracks that exhibit identical peak densities and therefore have 
constant expansion timescales. Our aim here is not a full parameter survey of conditions relevant to ^'^Ti 
synthesis, but rather an exploration of the sensitivity of ^^Ti synthesis to variations in crucial reaction rates 
affecting its production (and destruction) under conditions relevant to recent models and theory. Lacking 
detailed composition and velocity structure information our calculations cannot suggest firm predictions of 
^^Ti synthesis for comparison to observation. Rather this limited survey will take a more general approach 
that could serve to pla ce reasonable limits on the production of ^'*Ti with respect to 56,57,58j^j basis 
of nuclear systematics (IWooslev et al.lll99ll) . 



Our results will be given in terms of normalized production factors (NPF) P44, P57, and P58, that 
we define as the production factor for the given nucleus (the final mass fraction of the species in question 
divided by the mass fraction to which it decays in the Sun), normalized to the production factor for ^^Fe. 
Defined in this way NPF's for mass 56 (always made in these expansions predominantly as ^^Ni) are unity. 
Hence any NPF that is greater than unity will be overproduced with respect to solar iron, any NPF less 
than unit y will be underproduced. Since roug hly one-third to one-half of the iron in the Sun is attributable 
to SN II (jTimmes, Wooslev. fc Weaver] llQQsl ). NPF's for radioactive ".ssjs^j ^ ^^^^^^ of 2 - 3 would 

be acceptable (i.e. would not violate the solar ratio of ^''Fe/^^Fe). We consider these as upper bounds on 
the synthesis of ^^'^^Ni. They also serve to illustrate the range of electron fraction (Ye) allowed in such 
expansions. For referenc e, XC^'^Ca)fT, = 1.69 x 10-^ X(^*^Fe)0 = 1.26 x 10-^ X{^'^Fe)Q = 2.96 x 10-^ 
X(58Ni)o = 5.52 X 10-5 dLodderslbooi) . 



The decay of ^^'^^Co (with half-lives of 77.2 and 271.7 days, respectively) powers the light curves of 
SN II after the initial hydrogen recombination phase until late times when t he longer lived ^"''Ti and ^"^Co 
(with half-lives of 58.9 and 5.27 years, respectively) are expected to dominate (jTimmes et al.lll996[ ). Theory 
predicts that '^'^Ti is always produced in much smaller amounts than ^^'^''Ni, which we also infer from 
observations of Cassiopia A. This sets an upper bound on P44. Of course it could always be less, especially 
if SN II are not the only source of ^''Ti. 



4.2. 44rpj Nucleosynthesis in Select Supernova Models 



Before presenting our results we consider production o f ^^Ti and ^^■^'^ ■ ^^Ni f r om published m o dels o f 
SN II. The data a r e dra. wn fr om tables of stellar yields by Wooslev et al. ( 1995h . Rauscher et al. (2002), 
Limongi fc Chieflal(|2006h . and lYoung et al.l |2006l ) . We calculate mass fractions by normalizing the radioac- 
tive yields (in Mq) to the total mass of material ejected (initial mass — mass loss — remnant mass) and 
then form the normalized production factor. 

Figure [T2I shows the NPF's for each nucleus versus initial stellar mass for 11 < M/Mq < 25. We restrict 
our consideration to this mass range for two reasons: (1) for these solar metallicity models, the different 
treatments of mass loss are less pronounced for the lower masses, and (2) each survey explores a range of 
explosion energies for models with M/Mq > 30. Both can have a major impact on the ej ected mass and 
the m ass cut at higher masses. We do however show two results for the 25 Mq model of ^Rauscher et al" ^ 



(|2002[ ) where the explosion energy was increased to produce a model with double the yield of '^^Ni. We also 
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present results for a rang e of explosion energies in the 23 Mq model of I Young et al.l (|2006l ). The yields of 
Limongi fc ChiefR I ( 2006f) were derived assuming a constant '"^^Ni yield of 0.1 Mq. 



For ",58]vjj ^j^g results agree with the constraints imposed by the solar abundances, i.e. production factors 
between roughly 1 and 3 normalized to the production of ^^Fe in the Sun. The trends across the mass range 
shown, especially for ^^Ni, are quite similar, with the exception of the 20 Mq model of lRauscher et al.l ( 20021 ) 
that appears anomalous for reasons well described in that work. 

For ggg faij-iy remarkable consistency betwe en the three large surveys for the masses between 

18 and 25 Mq with NPF's ranging between 0.1 and 0.2. iLimongi fc ChiefR I (|2006f ) never produce it above 
0.1 for any mass, although its u niformity is in p art due to the constant ^^Ni yield ejected. They are however 
very similar to those of R aus cher et al.l (|2002i ). Of note is the larger value for P44 in a 23 model of 
Cassiopia A ( Young et al. 20061 ). The data reflect only their models with hydrogen envelopes removed to 



mimic a common-envelope evolution that also assumed a parameterized asymmetry of the explosion (those 
without asymetry ejected very little ^^Ni). The error bar reflects the variation of P44 due to their range 
of simulated explosion energies. A mixing algorithm was also included (absent in the other three surveys) 
whose role in the higher production of ^^^pj compared to the other surveys is unclear to us. 

Even at a factor of 2 agreement this uniformity may seem surprising for models that have very different 
prescriptions for many important physi cs elements like mass l oss (Woosley & Weaver had none), convection, 



opacities, and reaction rate libraries. IWoos 



Hauser-Feshbach reaction rate compilation of 



ev et al. ( 1995) was the only one to not use as a base the 



Rauscher fc ThielemannI |2000l ). 



More important still for "^^Ti production is the parameterization of the explosion and determination 
of the mass cut. All three of the large surveys use a paramet erization assuming motion of an inner zon e 
(i.e., a pisto n appr oach) as opposed to thermal energy input ( Thielemann. Nomoto. fc Hashimoto 19961 ). 
Young et al.l (|2006l ) vary their explosion energies through a parameterization of flux-limited diffusion for 
neutrinos that includes a "trapping radius" at which the neutrino opacity and flux is artificially adjusted. It 
should be noted that all these efforts initiate t he explosion a nd fol low the nuclear burning and subsequent 
ejecta in one-dimension with the exception of Young et al. ( 20061) who, after 10-100 s, map their results 
from one-dimension to three-dimensions to follow the mixing and determine the ultimate ejecta distribution. 
Unfortunately, their nucleosynthesis results are preliminary and do not include results for ^^'^^Ni, which 
should appear in a forthcoming work. 



4.3. 44rjij Nucleosynthesis Sensitivity Survey 

We now consider results for each of our simulated expansions defined in Table [H Figure [13] shows 
normalized production factors P44, P57, and P58 versus electron mole number Ye for adiabatic freeze outs 
fro m peak conditions defin ed for points CasA 1 — 3 in Table [S] All were drawn from a model for Cassiopia 
A ( Magkotsios et al. 20081) . In the figure, each central point represents a calculation that utilizes our rec- 



ommended ''^°Ca(a, 7)'^''Ti (production) rate for three choices of '*'*Ti(Q!,p)'*''V (destruction) rate. Solid line 
type and filled squares represent our "recommended" destruction rate, filled triangles represent its upper 
(dotted) and lower (dashed) bound. The error bars on each central point for all three surveys reflect the 
minimum and maximum deviations of P44 due to the six other choices of the '^^Ga,{a, 7)'*^Ti reaction rate that 
we considered. Tabulated nucleosynthesis results assuming our recommended semi-experimental production 
and destruction rates are given in Table |6l 
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The results for CasA points 1 and 2 are very similar. The NPF for ^^Ni and solar abundances suggest 
that the range of allowed electron fraction is 0.4980 < Ye < 0.500. As expected from nuclear systematics, in 
this range ^"'Ti is always less (by mass fraction) than ^''■^^Ni. It is also always underproduced with respect 
to solar iron (P44 < P^q = 1.0). An increase in P44 is seen at lower Ye, but this is due to a disparate drop in 
both 56Ni and ^-^Ti production (for point 2 between = 0.49 and 0.485, (X^^Ni) dropped by a factor of 46, 
while (X^^Ti) dropped by a factor of 3, respectively). For our "recommended" reaction rates the normalized 
production factors are very similar (~ 0.1 — 0.2 over the entire allowed range of Ye) to those seen in the SN 
II models with masses between 15 and 25 Mq (Figure fT2|). For the lower bound on our destruction rate, 
they are a factor of 1.7 higher. 

For point CasA-3 the results are very different. Now any choice of experimental rate for ^'^Ca(a, 7)^^Ti 
and either choice of our "recommended" or lower bound for the ^'*Ti(Q!,p)^''V destruction rate provides for 
nucleosynthesis that makes jj-^ proportion to solar iron over the allowed range of Ye- Integrated results 
from a SN model that experiences these conditions in a sizable fraction of its ejecta would be of interest. 
Adopting our upper bound for **^Ti(a,p)^''V still makes P44 between 0.2 and 0.8, or up to double the value 
seen in the two older stellar model surveys. For Ye less than 0.4980 the normalized production factors 
continue to chmb but are suggesting a strong overproduction with respect to solar iron (P57 > 3, Table O. 
For Y"e above 0.5, P44 is always small, even though the ^^Ni mass fraction remains high (^ 0.6). We also note 
that all three of these expansions were very a-rich (the final a-mass fractions across the Y"e range illustrated 
were 0.20, 0.23, and 0.32 for points CasA 1-3, respectively). 



The contour plot of Ti production in lMagkotsios et al.l ()2008[ ) suggests declining production for lower 



peak densities along the track of the CasA model. These translate to higher entropies (approaching 50 at 
Tgp = 4.2, pp — 5.0 X 10^ g cm~'^), but the temperature is dropping so fast that NSE (needed to make ^^Ni) 
will be increasingly hard to achieve. For these peak conditions, assuming Yg = 0.498 and our recommended 
principal production and destruction rates, the final a, ^^Ti, and 56,57,58]\jj mass fractions are 0.45, 1.8 x 10"'^, 
0.38, 0.03, and 0.08 respectively, leading to NPF's P44 = 3.53, P57 — 3.6, and P58 — 4.8. This again suggests 
a strong overproduction of the later two with respect to solar iron. 

For any given expansion the spread in values for P44 versus Yg is straight forward to understand. In 
general, the value of Yg tends to steer the net nuclear flows along pathways that are either proton-rich 
(Y"e > 0.5) or neutron-rich (Y"e < 0.5). Since ''^Ti and '^^Ni both have Y"e = 0.5, expansions of material that 
deviate far from this will reflect a drop in the abundance of each, alth ough for ^^Ni, which is made in NSE, 



the change will be less dramatic, especially for proton-rich conditions (jSeitenzahl et al.ll2008f ). For all three 
points (CasA 1-3) the principle production rate affecting ^'*Ti for all values ofY^ explored is ^°Ca(Q!, 7)^'*Ti. 
Differences in P44 across the range of Ye surveyed are therefore driven principally by the reactions affecting 
the production of ^"^Ca and the destruction of ^"^Ti. 

For 0.49 < Ye < 0.5 the dominant sequence of reactions producing '"^Ca are a-capture reactions on self- 
conjugate nuclei up to ■^^Ar followed by '^^ Ar{a, p)^^K{p, j)*'^Ca, operating principally between 2 < Tg < 3. 
This is also true for Ye > 0.5, but proton-induced (p, 7) and (j),a) flows provide additional pathways to 
build up '^^Ca at the expense of ^"^Ca. Proton-capture reactions also compete strongly with '**Ti(Q;,p)'*^V to 
accelerate the depletion of '^'^Ti over the same temperature range. For Ye < 0.49 proton-induced reactions 
that produce ^"Ca and destroy ^^Ti are suppressed due to the lower free proton abundance, the ones affecting 
production win out. 

As mentioned above, the error bars on each point represent the variation in P44 due to the many choices 
of '*°Ca(Q!, 7)^'*Ti reaction rate we considered. For point CasA-3 this amounted to roughly 20% deviations 
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from our "recommended result", most often bracketed by our x2 upper and lower error bar. Considering 
only the rates from other efforts, the variation would have been +20% and -10% from our recommended 
value. The variation is asymetric due to the relative ratios of the various production rates versus ours which 
was used to anchor the central point (see inset of Figure [S]) ■ If instead of the factor of 2 errors we assumed 
for the supplemental theory cross section we had used the errors suggested by our off-line counting data 
(Figure |4]), the variation would have been of order 2%. 

The choice of '*^Ti(a,p)^''V destruction rate exhibits a larger sensitivity (^ 70%) to overall ''^Ti synthesis 
than the entire spread in the experimental '^^Ca,{a, 7)^*Ti production rates. Taken together, the potential 
uncertainty over the allowed range in ^''Ti synthesis encompassed by the possible choices of production 
and destruction rate is roughly a factor of 3. 

Interestingly, t he dot-short-dash lines in Figu re fTSl show ''^Ti synthesis considering the use of the orig- 
inal published fits ( Rauscher fc Thielemann 2000) for both '*°Ca(a, 7)^'*Ti and '*'*Ti(a,p)*^V, in very close 



agreement (^ 4%) with what resulted using our "recommended" rates. The reason is that for 1 < Tg < 4, 
the original fit to the production rate (Figure IH]) made on average (more above Tg ~ 2.4, less below) very 
nearly the same amount of ''^Ti as our "recommended rate", while the destruction rate (Figure [TU]) was 
nearly identical for Tg > 3. Our higher production rate below Tg = 2.4 accounted for the uniformly higher 
production. Again, this is the same level of variation from our recommended rate that would have occurred 
if we had considered only our offline counting error (Table |4]). This would also be the result obtained using 
the default BDAT reaction rate library provided with the TORCH reaction network code (Timmes 2010, 
private communication) , since i t utilizes the same (ori ginal) fits to the relevant H auser-Feshbach theory rates 
( Rauscher fc Thielemann 2000 ') and 3a reaction rate ( Caughlan fc Fowler 19981 ). 



Figure [14] shows our results for points 2DExpl 4 — 6 (Table [6]) drawn from conditions in a rotating 
two-dimensional SN explosion model. These have the same peak temperature albeit at hi gher peak density, 



as the points drawn from the CasA model, and sample the same regions as depicted in iMagkotsios et al 



(|2008[ ). The results for P57, and P58 are very similar to those from the model for CasA (i.e. production 
compared to solar iron), but P44 is much smaller for all scenarios. It is similar however to that seen in the 
massive star surveys (P44 ^ 0.1 — 0.2, Figure [T2|). The final a-particle mass fractions across the Ye range 
illustrated were roughly 0.10, 0.07, and 0.06 for points 4-6, respectively (Table |6|). 

According to^M agkotsios et"aL ( 20081 ) points 2DMIID 7 — 8 represent conditions in a model for a rotating 



two-dimensional MHD star. Point 7 is located in the "Chasm" of ^^Ti production, while point 8 is in the 
region defined by them as incomplete silicon burning. Both experience much higher peak densities (at 
identical peak temperatures) than points from the previous two models, and consequently they have very 
low entropies and very short hydrodynamic time scales. Results assuming our "recommended" production 
and destruction rates are shown in Table [6l As before P57 and P58 are very similar to our previous results 
reflecting a near constant production of ^'^'^^Ni with respect to solar iron (for a given Yg) over the entire 
range of expansion timescales surveyed. However, P44 was virtually non- existent, with (X^^Ti) < 10^'^ 
(often much less) for any Yg in either expansion. Consid ering point 7, whose peak conditions reside in the 



"Chasm" of ■*'*Ti production noted in iMagkotsios et al. ([2OO8), our "recommended" reaction rates would 



slightly lessen (by 4%) the "depth" of the chasm compared to the reaction rates they used. 

The explanation for the drop in P44 stems from the lower entropy (higher peak density) for each point in 
these rotating two-dimensional explosion model expansions versus those from the model for CasA (see Table 
[6|). Similarly, for all three CasA points, the lower entropy expansions translated into a lower P44 (Figure 
[13]). At lower entropy material merges very quickly into the iron group at the expense of the a-particle 
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abundance which at late times was so low (~ 1%) in the two-dimensional explosion model expansions that 
few were available to make ^°Ca at appreciable levels, with 44^j being consequently even lower. 

From the track of the rotating two-dimensional explosion model we chose conditions (points 2DExpl 
4 — 6) that had identical peak densities, and hence a constant expansion timescale (0.14 s). For this factor 
of 2.6 overall change in entropy and for the specific case of Ye = 0.4980, the final ^^Ni mass fractions were 
0.74, 0.77, and 0.77 for points 2DExpl 4 — 6, respectively. For the CasA model, with a factor of 1.5 range in 
entropy (that included a factor of 2 range of dynamic time-scale) the ^^Ni mass fractions were 0.64, 0.62, and 
0.52 for points CasA 1-3, respectively. The salient point for '^''Ti synthesis again is that the higher entropy 
expansions made less ^^Ni and had more a-particles, especially at late times, that enabled the production 
of ""^Ca and '''^Ti (Table[6]). Since P44 is a ratio of Ti to Fe, it grows as a function of increasing entropy. 

We illustrate all these points graphically in Figure [15] where we show the evolution of select mass 
fractions versus temperature for the three expansions in our survey with identical peak temperatures (points 
2, 5, and 8) that all started with an initial composition with — 0.4980. For each expansion both ■*°Ca and 
^"^Ti achieve an NSE abundance not very different from their final freeze out values, but both are effectively 
destroyed for 4 < Tg < 5. Thereafter they are reassembled at a rate determined largely by the a-particle 
abundance, which is affected by the peak density (entropy). It is clear that a more robust a-rich freeze out 
is more conducive to ^^Ti synthesis and that reaction rates clearly matter. 

We also note for all conditions shown in Figure [15] the freeze-out of the ^*Ti abundance near Tg — 2. 
The ''^V abundance does increase as the temperature declines, but the '*^Ti abundance changes by less than a 
few percent (the increase is due predominantly to the decay of radioactive ''^Cr). At this Ye and for Tg > 2.0 
the dominant destruction rate for ''^Ti is ''^Ti(Q!,p)'*^V, below 1.0 it switches to '*^Ti(p, 7)*^V. Terminating 
our simulations at any point below Tg — 2.0 would have a negligible effect on P44. 



4-. 3.1. Effect of the Expansion Timescale 

''■^Ti synthesis is also affected by the expansion timescale, thd — ^"^^xl \Jlh>- Figure [16] shows the range 
of nuclei produced in an a-rich freeze-out from peak conditions given for point CasA-3 in Table [5] with an 
initial composition of Ye — 0.4980 and our default scaling of % — 1. The nucleosynthesis is presented in 
terms of "traditional" production factors (so far our results have been expressed as ratios of production 
factors to that of ^^Fe). In the figure isotopes of a given element are connected by solid lines, those produced 
as radioactive progenitors are surrounded by a diamond, the most abundant isotope of a given element is 
denoted by a star. The dotted lines indicate a factor of 2 above and below the dashed line centered on ^^Fe 
(made as ^^Ni, ''^Ca is made as ^''Ti). The dominant species are in the iron group. This is a fairly typical 
result for any of the CasA or two-dimensional explosion model expansions within the allowed range of Ye . 

Figure [17] explores the effect of increasing the hydrodynamic time-scale by a factor of x = 5. Shown 
is a straight forward ratio of " traditional" production factors for the CasA-3 expansion assuming both the 
default (x = 1) and extended (x — 5) scaling on thd- The species in the iron group are not much affected, 
56,58]\jj ^j.g siig]-^tly higher than unity, all others show lower production factors. The reason is that for a 
longer expansion timescale, the material experiences higher temperatures for a longer duration. Those made 
dominantly in NSE (the iron group seen here) are slightly altered due to a-particle reactions on them, which 
are ultimately provided by photo-disintegration and operation of the 3a-reaction. The effect is to reduce 
the overall fraction of a-particles available at late times for the re-assembly of species lighter than the iron 
group. 
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Figure [18] shows for the specific CasA-3 example above the temperature evolution over the first 4 s 
assuming the two scalings. Over that time the temperature in the default (x = 1, dotted line) scenario 
declined from Tg = 4.7 - 0.25 while the longer expansion (x = 5, dashed line) only declined to Tg — 2.6. 
Also shown is the ratio of the a-particle mass fraction (long/default). Recall that both ^"Ca and ''^Ti begin 
to re-assemble below Tg = 4 (Figure [15]) . By the time the longer expansion scenario reaches Tg = 4, it only 
has 65% of the a-particle mass fraction that the default expansion had, by Tg = 3 it is down to 52%. With 
only half the a-particles to work with, the expansion with the higher scaling only reached P44 — 0.5 ±0.1, a 
factor of 2.6 less than the default expansion. Figure [T9l shows the effect over all scenarios considered in the 
CasA model expansions. The average P44 for points 4-6 (the two-dimensional explosion model) were 0.1, 
0.05, and 0.08, respectively, over the allowed range of Y^- For points 7 and 8 they are essentially zero. 



5. Conclusions 

We have considered the sensitivity to 44^j production in expansions that approximate freeze outs from 
NSE due to variations in the principle production and destruction reactions ^°Ca(Q;, ^Y^'Y'y and ''**Ti(a,p)''''V, 
and contrast them to experimental and theory reaction rate developments over the past 20 years. Experimen- 
tally we have also measured a thick-target yield for the ^'^Ca(Q;,7)^'*Ti reaction. In-beam 7-ray spectroscopy 
was used to determine the yield of the 1083 keV prompt 7-ray from ''^Ti at = 4.13, 4.54, and 5.36 
MeV. In order to correct for those transitions which bypass the 1083 keV transition, the Monte Carlo code 
DICEBOX was used to estimate a correction of 20% to the in-beam thick target yield. An off-line activation 
measurement using the target from the = 5.36 MeV irradiation showed good agreement with the in-beam 
measurement. 

We then derived a thermonuclear reaction rate by normalizing the NON-SMOKER cross section (down 
by a factor of 1.71) to agree with our measured off-line thick target yield. We derived an error bar for 
our recommended rate whose magnitude (± x 2) was dominated by the theoretical cross section error. We 
als o carry out a similar re-evaluation of the ''^Ti(a,p)''''V reaction rate whose cross section was measured 



bv ISonzogni et al.l (|2000l ). For both reactions, we conclude that the experimental data were far above the 
Gamow window, and suggest that further measurements be attempted. 

We then carried out a sensitivity survey of ^^Ti nucleosynthesis in adiabatic expansions from eight 
peak temperature and d ensity combinations drawn from conditions in three recent stellar explosion models 
( Magkotsios et al. 20081 ). For each expansion we survey a range of initial compositions (0.505 > Ye > 0.485). 



We also vary the principle production and destruction rates affecting '^^Ti in these expansions (eight choices 
of production rate, four for destruction). Our results show ^^^j produced in proportion to solar iron for only 
one expansion drawn from a model for Cassiopia A, even though the final mass fractions for 44^j ^nost 
of the expansions are typically of order 10^"*. With one exception, the other expansions are consistent with 
^"^Ti production seen in previous surveys of one-dimensional stellar evolution and with constraints imposed 
by solar abundances. Our results suggest that a strong a-rich freeze out {X{a)f ^ 0.2 — 0.3) is highly 
conducive to '^'^Ti synthesis. 

With respect to reaction rate sensitivity, our experimental results suggest a recomme nded '^^Caja. 7)'^^Ti 



react i on rate that is sma ller than those predicted by the most recent experimental efforts (jVockenhuber et al 



2007; Nassar et al. 2006h . but with a fairly large error bar (± x 2) that would in fact encompass the former 



one. Nucleosynthesis models using our recommended rate would suggest less **Ti tha n these two recent 



experiments, although it would be higher than that suggested by an earlier measured rate ([Cooperman et al 
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1977 ) and the current theory rate ( Rauscher fc Thielemannlboool) . The total range of sensitivity is a factor of 
1.5 (36%) when considering all of the production rates available. We also find that the uncertainties associated 
with the dominant destruction rate, '*'*Ti(a,p)'*''V, have roughly double the impact on ''^Ti synthesis (x3.2, 
or 70%) than that exhibited by the entire range of production rates. W e suggest the use of our re-evaluation 
of the only available experimental reaction rate ( Sonzogni et al. 200Cll) and strongly suggest consideration of 
its attendant larger uncertainty (± x 3). This is a sli ghtly larger spread in '^'^Ti sensitivity t han observed 
in two recent surveys of massive star nucleosynthesis (jRauscher et al.ll2002t iLimongi fc Chief fi 2006) that 
used essentially the same nuclear data, suggesting that current uncertainties in reaction rates could lead to 
as large an uncertainty in 44^j synthesis as that produced by different treatments of stellar physics. 



Since the seminal work of IWooslev et al.l (jl973[ ). several new and novel theories of SN nucleosynthesis 
have featured regimes withi n the neutrino wind where the a— rich freeze out plays an importa nt role, including 



scenarios for the r-proces s (jWooslev et al. Ill994i iHoffman et al.lll997l ). and the j/p-process (jFrohlich et. al. 
2006HPruet et. al. hooj] . In each a combination of high entropy (S'rad ^ 50) and short expansion time- 



scale are required to produce the unique signatures of each process (the r-abundances, and light p-nuclei, 
respectively). However, as a consequence of extreme neutrino irradiation, the composition is forced away 
from neutron-proton equality and ultimately high entropy material enters regions above the iron group where 
local effects due to rapid changes in particle separation energies have a strong influence on the net nuclear 
flows. In these works no ^^^j jg reported. 

F urther out in the SN electa theory does no t show ^''^Ti enhancement due to the z^-process in massive 



stars |Wo oslev et al. 



scenarios (Hoffma n et al 



1990; 

T 



2008 



Wooslev et a l.lll995l). nor in recent low-mass (electron capture) core-collapse 
Wanajo e t al. 2009) ■ The later has been shown to effectively synthesize the 
long elusive a-rich freeze-out candidate ®'*Zn, but not ^*Ti. Our limited survey suggests that '^^Ti synthesis 
requires modest entropy (iSrad ^ 35) and expansion timescales (^ thd = 0.45 s) over a fairly narrow range 
of Ye (0.5 >Ye> 0.4980). 

To date surveys of massive star evolution and nucleosynthesis have typically underprodu ced species 
whose production is attributed to the a-rich freeze-out, in particular ^''Ca (made as ^"^Ti) and ^*Zn (jWooslev et al 



1995t iLimongi fc Chiefln 1120061 ). This has largely been ascribed to variations in treatments of stellar physics. 



most notably the parameterization of the explosion in one-dimensional models. Incorporating stellar yields 
from these surveys into models for galactic chemical ev olution indicate the degree of underproduction, roughly 
a factor of a 2 — 3 (jTimmes. Wooslev. fc Weaveiill995[ ). Recent work has lead to models for the progenitor of 
Cassiopia A that suggest increased production, but the underlying physics is still quite uncertain and in need 
of improvement ("Youn g et all2006l) . Future attention may focus on issues related to asymmetrical explosions 
where a noted increase in " ^"^Ti production compared to models with i mposed spherical symm e try has been 



1997 



Nagataki et al. 19981 : Hwang fc Laming 20031 : Young et al 



sugge sted for many years (jNagataki et al 
20061 ). The community is on the threshold of three-dimensional calculations that should provide valuable 



insight into the core-collapse mechanism, including physics such as rotation and magnetic fields which will 
likely enforce an asymmetrical result. Such results should help guide future parameterizations of the explo- 
sion in our one-dimensional models which will continue to carry the burden in future surveys of massive star 
evolution and nucleosynthesis. We believe we have addressed the uncertainty in two key nuclear reaction 
rates affecting "'^Ti synthesis, but ultimately, if Type II SNe are the dominant site of '*^Ti production, future 
models will have to include a larger fraction of their ejecta that experience an a-rich freeze out than they 
have in the past. Another solution would be an additional s ource of ^"^Ti, such as rare SNe of type la 
IWooslev et al.lll986t Iwooslev Ill997t ) or lb |Perets et al. II2OI0I ). 
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Table 6. Nucleosynthesis Survey 



Ye 


0.5050 


0.5000 


0.4995 


0.4990 


0.4985 


0.4980 


0.4965 


0.4950 


0.4900 


0.4850 


v 


-0.0100 


0.0000 


0.0010 


0.0020 


0.0030 


0.0040 


0.0070 


0.0100 


0.0200 


0.0300 




Point 1: Model for Cassiopia A, Tg^ 


= 6.5, pp = 


4.0 X 10^ Srad = 22. 


8, thd = 


0.22 s 




X{a) 


2.06(-l) 


2.08(-l) 


2.07(-l) 


2.06(-l) 


2.04(-l) 


2.03(-l) 


1.98(-1) 


1.94(-1) 


1.79(-1) 


1.64(-1) 


X(40Ca) 


1.66(-4) 


8.27(-4) 


7.87(-4) 


7.46(-4) 


7.15(-4) 


6.89(-4) 


6.26(-4) 


5.72(-4) 


4.25(-4) 


1.00(-4) 


X(44Ti) 


6.17(-5) 


3.00(-4) 


3.49(-4) 


3.42(-4) 


3.34(-4) 


3.27(-4) 


3.05(-4) 


2.85(-4) 


2.22(-4) 


7.39(-5) 


X(56Ni) 


7.51(-1) 


7.35(-l) 


7.17(-1) 


6.92(-l) 


6.66(-l) 


6.39(-l) 


5.58(-l) 


4.77(-l) 


2.11(-1) 


4.58(-3) 


X('"Ni) 


2.09(-3) 


1.22(-2) 


2.89(-2) 


3.09(-2) 


3.27(-2) 




3.78(-2) 


3.98(-2) 


3.48(-2) 


4.52(-3) 


X(5SNi) 


2.25(-5) 


1.57(-3) 


9.70(-3) 


3.50(-2) 


6.04(-2) 


8 5q('-2'l 


1.63(-1) 


2.41(-1) 


5.08(-l) 


6.97(-l) 


P44 


0.062 


0.305 


0.364 


0.369 


0.374 


0.382 


0.409 


0.446 


0.786 


12.027 


P57 


0.131 


0.707 


1.717 


1.891 


2.079 


2.283 


2.889 


3.536 


7.024 


42.192 


P58 


0.004 


0.221 


0.309 


1.153 


2.060 


3.071 


6.659 


11.504 


54.762 


3452.055 




Point 2: Model for Cassiopia A, Tgp 


= 5.5, pp = 


2.0 X 10^, 


Srad — 27.7, ThD = 


0.32 s 




X{a) 


2.22(-l) 


2.25(-l) 


2.25(-l) 


2.23(-l) 


2.22(-l) 


2.21(-1) 


2.16(-1) 


2.12(-1) 


1.98(-1) 


1.84(-1) 


X(40Ca) 


2.47(-4) 


9.50(-4) 


9.06(-4) 


8.61(-4) 


8.27(-4) 


7.99(-4) 


7.30(-4) 


6.73(-4) 


5.15(-4) 


1.04(-4) 


X(44Ti) 


8.66(-5) 


3.59(-4) 


4.14(-4) 


4.06(-4) 


3.98(-4) 


3.90(-4) 


3.68(-4) 


3.46(-4) 


2.80(-4) 


8.13(-5) 


X(56Ni) 


7.35(-l) 


7.17(-1) 


6.99(-l) 


6.75(-l) 


6.48(-l) 


6.22(-l) 


5.40(-l) 


4.59(-l) 


1.92(-1) 


3.06(-3) 


X("Ni) 


2.18(-3) 


1.22(-2) 


2.82(-2) 


3.00(-2) 


3.17(-2) 


3.32(-2) 


3.65(-2) 


3.82(-2) 


3.25(-2) 


3.17(-3) 




2.53(-5) 


2.03(-3) 


1.00(-2) 


3.54(-2) 


6.09(-2) 


8.64(-2) 


1.64(-1) 


2.42(-l) 


5.09(-l) 


6.57(-l) 


P44 


0.088 


0.374 


0.441 


0.450 


0.458 


0.468 


0.508 


0.562 


1.085 


19.836 


P57 


0.132 


0.721 


1.715 


1.884 


2.078 


2.267 


2.867 


3.534 


7.190 


44.672 


P58 


0.004 


0.216 


0.328 


1.198 


2.136 


3.178 


6.900 


12.000 


60.261 


4877.049 



Point 3: Model for Cassiopia A, Tgp = 4.7, pp = 1.0 x 10^, Srad =■ 34.5, thd = 0.45 s 





3.13(-1) 


3.17(-1) 


3.16(-1) 


3.15(-1) 


3.14(-1) 


3.13(-1) 


3.09(-l) 


3.06(-l) 


2.96(-l) 


2.86(-l) 


X(40Ca) 


7.70(-4) 


2.10(-3) 


2.04(-3) 


1.94(-3) 


1.87(-3) 


1.82(-3) 


1.69(-3) 


1.59(-3) 


1.31(-3) 


1.85(-4) 


X(44Ti) 


2.05(-4) 


8.57(-4) 


9.81(-4) 


9.69(-4) 


9.57(-4) 


9.45(-4) 


9.10(-4) 


8.77(-4) 


7.74(-4) 


1.44(-4) 


X(56Ni) 


6.40(-l) 


6.17(-1) 


5.99(-l) 


5.75(-l) 


5.49(-l) 


5.22(-l) 


4.41(-1) 


3.60(-l) 


9.43(-2) 


2.56(-3) 


X(57Ni) 


1.83(-3) 


1.27(-2) 


2.81(-2) 


2.99(-2) 


3.15(-2) 


3.29(-2) 


3.55(-2) 


3.62(-2) 


2.39(-2) 


1.30(-3) 


X(58Ni) 


2.21(-5) 


2.56(-3) 


8.26(-3) 


3.29(-2) 


5.77(-2) 


8.26(-2) 


1.58(-1) 


2.34(-l) 


4.94(-l) 


4.31(-1) 


P44 


0.238 


1.037 


1.223 


1.258 


1.300 


1.349 


1.543 


1.818 


6.120 


41.053 


P57 


0.123 


0.878 


1.994 


2.210 


2.454 


2.675 


3.429 


4.266 


10.787 


22.010 


P58 


0.003 


0.243 


0.315 


1.304 


2.408 


3.614 


8.171 


14.825 


119.333 


3732.058 


Point 4: Rotating Two-dimensional Explosion Model, Tgp 


= 6.5, Pp 


= 1.0 x 10^, 


Srad = 9.14, THD = 


0.14 s 


X[a) 


1.03(-1) 


1.06(-1) 


1.05(-1) 


1.04(-1) 


1.03(-1) 


1.02(-1) 


9.77(-2) 


9.39(-2) 


8.14(-2) 


6.95(-2) 


X(40Ca) 


2.00(-5) 


4.54(-4) 


4.32(-4) 


4.03(-4) 


3.81(-4) 


3.62(-4) 


3.15(-4) 


2.75(-4) 


1.74(-4) 


8.54(-5) 


X(44Ti) 


3.11(-5) 


1.34(-4) 


1.90(-4) 


1.83(-4) 


1.77(-4) 


1.70 (-4) 


1.51(-4) 


1.34(-4) 


8.76(-5) 


4.87(-5) 


X(56Ni) 


8.50(-l) 


8.40(-l) 


8.21(-1) 


7.94(-l) 


7.67(-l) 


7.39(-l) 


6.56(-l) 


5.74(-l) 


3.04(-l) 


4.68(-2) 


X(57Ni) 


4.15(-3) 


1.13(-2) 


2.79(-2) 


3.07(-2) 


3.32(-2) 


3.55(-2) 


4.07(-2) 


4.40(-2) 


4.30(-2) 


1.98(-2) 


X(58Ni) 


7.02(-5) 


9.96(-4) 


1.14(-2) 


3.58(-2) 


6.04(-2) 


8.51(-2) 


1.60(-1) 


2.37(-l) 


5.03(-l) 


7.79(-l) 


P44 


0.029 


0.119 


0.173 


0.173 


0.172 


0.172 


0.172 


0.175 


0.216 


0.777 


P57 


0.267 


0.571 


1.445 


1.648 


1.839 


2.044 


2.649 


3.268 


6.017 


18.011 


P58 


0.021 


0.183 


0.317 


1.029 


1.790 


2.624 


5.566 


9.408 


37.801 


379.032 



^ All results assume our "recommended" rates for ■*°Ca(Q:, 7)4'*Ti and '*4Ti(a,p)4^V. 
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Table 6. Nucleosynthesis Survey^ (continued) 

Ye 0.5050 0.5000 0.4995 0.4990 0.4985 0.4980 0.4965 0.4950 0.4900 0.4850 

T] -0.0100 0.0000 0.0010 0.0020 0.0030 0.0040 0.0070 0.0100 0.0200 0.0300 

Point 5: Rotating Two-dimensional Explosion Model, Tgp = 5.5, pp = 1.0 x 10^, Srad = 5.54, thd = 0.14 s 

X{a) 6.85(-2) 7.23(-2) 7.17(-2) 7.07(-2) 6.96(-2) 6.85(-2) 6.53(-2) 6.21(-2) 5.23(-2) 4.36(-2) 

X(40Ca) 1.13(-5) 3.68(-4) 3.56(-4) 3.29(-4) 3.08(-4) 2.90(-4) 2.46(-4) 2.10(-4) 1.23(-4) 6.11(-5) 

X(44Ti) 1.61(-5) 8.16(-5) 1.62(-4) 1.54(-4) 1.47(-4) 1.41(-4) 1.22(-4) 1.05(-4) 6.34(-5) 3.43(-5) 

xl^^Wi) 8.77(-l) 8.74(-l) 8.54(-l) 8.27(-l) 7.99(-l) 7.71(-1) 6.87(-l) 6.04(-l) 3.31(-1) 6.85(-2) 

X("Ni) 7.18(-3) 9.75(-3) 2.59(-2) 2.92(-2) 3.22(-2) 3.48(-2) 4.08(-2) 4.45(-2) 4.47(-2) 2.40(-2) 

Xf^Ni) 1.86(-4) 8.99(-4) l-18(-2) 3.54(-2) 5.92(-2) 8.33(-2) 1.57(-1) 2.32(-l) 4.95(-l) 7.73(-l) 

P44 0.015 0.070 0.142 0.139 0.137 0.136 0.132 0.130 0.143 0.373 

P57 0.451 0.474 1.291 1.502 1.717 1.928 2.527 3.125 5.741 14.926 

P58 0.056 0.157 0.316 0.977 1.685 2.467 5.201 8.771 34.144 257.353 

Point 6: Rotating Two-dimensional Explosion Model, Tgp = 4.7, pp = 1.0 x lO'^, Srad = 3.45, thd = 0.14 s 

X{a) 5.67(-2) 6.07(-2) 6.04(-2) 5.97(-2) 5.90(-2) 5.84(-2) 5.65(-2) 5.48(-2) 4.99(-2) 4.58(-2) 

X(40Ca) 1.36(-5) 4.91(-4) 5.04(-4) 4.67(-4) 4.42(-4) 4.20(-4) 3.71(-4) 3.32(-4) 2.36(-4) 1.46(-4) 

X(44Ti) 1.19(-5) 6.90(-5) 2.42(-4) 2.33(-4) 2.24(-4) 2.17(-4) 1.96(-4) 1.79(-4) l-34(-4) 9.15(-5) 

X(56Ni) 8.68(-l) 8.76(-l) 8.56(-l) 8.28(-l) 7.99(-l) 7.71(-1) 6.86(-l) 6.02(-l) 3.26(-l) 6.19(-2) 

X("Ni) 1.14(-2) 8.83(-3) 2.47(-2) 2.85(-2) 3.18(-2) 3.46(-2) 4.09(-2) 4.46(-2) 4.42(-2) 2.28(-2) 

Xlc'^m) 4.46(-4) 8.06(-4) 1.05(-2) 3.22(-2) 5.43(-2) 7.66(-2) 1.45(-1) 2.15(-1) 4.59(-l) 7.14(-1) 

P44 0.011 0.059 0.210 0.210 0.209 0.209 0.213 0.222 0.306 1.102 

P57 0.720 0.428 1.229 1.464 1.685 1.912 2.532 3.159 5.753 15.650 

P58 0.139 0.143 0.281 0.888 1.550 2.271 4.826 8.159 32.124 262.195 

Point 7: Rotating Two-dimensional MHD Model (CHASM region), Tgp = 6.5, pp = 1.0 x 10^, Siad = 0.91, thd = 0.04 s 

X{a) 1.01(-2) 1.31(-2) 1.20(-2) 1.08(-2) 9.65(-3) 8.51(-3) 5.25(-3) 2.26(-3) 1.81(-6) 8.95(-7) 

X(40Ca) 7.43(-8) 3.61(-5) 2.64(-5) 1.94(-5) 1.41(-5) 9.95(-6) 2.68(-6) 2.88(-7) 4.06(-8) 2.35(-7) 

X(44Ti) 1.22(-7) 1.15(-5) l.ll(-5) 8.25(-6) 5.97(-6) 4.19(-6) 1.08(-6) 1.06(-7) 1.98(-9) 6.31(-9) 

X(56Ni) 8.67(-l) 9.40(-l) 9.14(-1) 8.87(-l) 8.59(-l) 8.32(-l) 7.52(-l) 6.71(-1) 3.93(-l) 1.28(-1) 

X("Ni) 1.06(-2) 1.34(-2) 2.30(-2) 2.84(-2) 3.30(-2) 3.69(-2) 4.57(-2) 5.15(-2) 4.88(-2) 3.07(-2) 

X(58Ni) 7.53(-4) 1.50(-4) 1.92(-2) 4.21(-2) 6.61(-2) 9.11(-2) 1.72(-1) 2.60(-l) 4.30(-l) 5.61(-1) 

P44 0.000 0.009 0.009 0.007 0.005 0.004 0.001 0.000 0.000 0.000 

P57 1.327 0.607 1.069 1.365 1.642 1.891 2.596 3.265 5.272 10.388 

P58 1.244 0.064 0.478 1.082 1.760 2.496 5.209 8.837 24.888 99.029 

Point 8: Rotating Two-dimensional MHD Model, Tgp = 5.5, pp = 1.0 x 10^, Srad = 0.55, thd = 0.04 s 

X{a) 9.27(-3) 1.13(-2) 1.07(-2) 9.89(-3) 9.16(-3) 8.47(-3) 6.57(-3) 4.86(-3) 4.69(-4) 2.45(-6) 

X(40Ca) 1.39(-7) 5.28(-5) 4.50(-5) 3.56(-5) 2.84(-5) 2.27(-5) 1.12(-5) 4.92(-6) 1.59(-8) 4.70(-10) 

X(44Ti) 2.14(-7) 6.95(-6) 2.06(-5) 1.65(-5) 1.33(-5) 1.07(-5) 5.20(-6) 2.24(-6) 5.93(-9) 4.65(-ll) 

X(56Ni) 7.86(-l) 9.25(-l) 9.01(-1) 8.74(-l) 8.47(-l) 8.20(-l) 7.39(-l) 6.59(-l) 3.89(-l) 1.18(-1) 

X(57Ni) 1.07(-2) 8.79(-3) 2.10(-2) 2.67(-2) 3.14(-2) 3.52(-2) 4.39(-2) 4.95(-2) 5.37(-2) 3.37(-2) 

X(58j^i) 1.40(-3) 2.13(-4) l-47(-2) 3.50(-2) 5.64(-2) 7.89(-2) 1.52(-1) 2.32(-l) 5.45(-l) 7.60(-l) 

P44 0.000 0.006 0.017 0.014 0.012 0.010 0.005 0.003 0.000 0.000 

P57 1.505 0.404 0.992 1.298 1.575 1.828 2.521 3.193 5.890 12.301 

P58 2.660 0.086 0.373 0.914 1.516 2.197 4.685 8.031 31.974 146.341 
^ All results assume our "recommended" rates for ^°Ca(Q;, 7)'^^Ti and ^^Ti(a,p)^^V. 
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Fig. 1. — Partial decay schemes of Ti and its daughter "'^Sc. 




2. — A Schematic of the target chamber used in the experiment. 
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Fig. 3. — Partial HPGc 7-ray spectra at Ea = 5.36 for the detector at 99° with a simultaneous fit to the 
1039 keV ^"Ge and 1083 keV '''^Ti 7-rays. 
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Fig. 4. — Measured thick target yields compared to the calculated NON-SMOKER thick target yield. The 
online measurement is the measured yield for the 1083 keV 7-ray corrected upward by 20% (see the text). 
The Ea = 5.36 MeV offline counting data point and its error is plotted as measured. The NON-SMOKER 
cross section required a downward scaling of 1.71 to match the off-line data point. The dashed lines indicate 
the normalization needed to match the offline error bars. 
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Fig. 5. — 7-ray spectra observed in a two week low background count of the activated target bombarded at 
= 5.36 MeV. 
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Fig. 6.— The Effective Gamow windows for '^°Ga{a,-f)'^*Ti and '^*Ti{a,p)^'^Y at Tg = 2 and Tg = 5. Shown 
are the true Gamow window as given by the real integrand of the reaction rate formula (solid line type, see 
Eq. U) , the product of the high-energy tail of the Maxwell-Boltzmann distribution with the simple Coulomb 
barrier penetration factor (denoted as "exp(-E/kT-b/sqrt(E)), dot-dash line type)" and the usual Gaussian 
approximation to the latter ( "Gauss approx" , dotted hue type). Stars deno te the energy values at which 
experimental data were measured in our effort and from lSonzogni et al.l (|2000[ ). The Gaussian approximation 
is widely used to estimate the Gamow window but predicts it at too high an energy, especially for the capture 
reaction. For purposes of illustration, all curves were renormalized to the arbitrary unit scale shown, for 
Tg — 5 all three quantities peak at 10, for Tg = 2 they peak at 6. 
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Fig. 7. — Comparison of the experimental energy range for ^"Ca(a.7)'*^Ti with the true Gamow window 
given by the actual maximal contributions to the reaction rate integral. The experimental data do not 
contribute significantly to the rate in the temperature range relevant to ''^Ti synthesis (2.0 < Tg < 4.0). 
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Fig. 9. — Fits to reaction rates for ^°Ca(Q;,7)'*^Ti considered in this study. Solid and dashed hne types denote 
experimental rates from various authors, dotted and dot-dash line types are from various theory efforts. The 
tabulated NON-SMOKER theory rate (jRauscher fc ThielemannI 120011 ) is denoted by crosses. Our recom- 
mended rate (and its error, ± x 2) are indicated as solid (lighter) lines. The inset ill ustrates the ratio of each 
reaction rate to the NON-SMOKER rate, including the (dot-short-dash) original fit (jRauscher fc Thielemann 



20001). our attempt (dot ted) to refit the tabulated NON-SMOKER rate, and the theory rate (dot-long-dash) 



of 



Wooslev et alf|l978l) 
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Fig. 10. — Fits to reaction rates for '^**Ti(a,p)'*''V considered in this study. Solid li ne type denotes our re - 
evaluation and its likely error (± x 3) of the lone ex perimental rate (dashed-line typ e) of lSonzogni et al.l (|200(3l ) , 
while the tabulated NON-SMOKER theory rate (jRauscher fc Thielemannll2001l ) is denoted by crosses. The 
inset illustr ates the ratio of each reaction rate to the NON-SMOKER rate, including the (dot-short-dash) 
original fit ( Rauscher fc Thielemannll2000l) , our attempt ( dotted ) to refit the tabulated NON-SMOKER rate, 
and the theory rate (dot-long-dash) from I Wooslev et aL ( IQTSi ). 
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Fig. 11. — Re-evaluated ^^T i(a,p)^''V rate, based on data by lSonzogni et al.l (|2000l ): Ratio of the new rate 
to the NON-SMOKER rate (jRauscher fc ThielenianrJl200lh . The estimated range of the rate is derived from 
the theoretical and experimental errors (see the text for details). The geometric mean of the upper and lower 
limit is used as a recommended rate in this work. The impact of the Sonzogni et al. ( 2000l ) data manifests 
itself in the deviation of the geometric mean from unity. Within errors, the new rate is compatible with the 
NON-SMOKER rate. 
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Fig. 12.— Normalized production fa ctors for ^"^Ti, ^^Ni and ^^Ni versus initial stellar mass derived from 
nucleosynthetic yi elds of SN models ( Wooslev et al. IQQSt Rauscher et al. 2002 : Limongi fc ChiefR 2006t 



Young et al.ll2006l ) 



- 39 - 




Fig. 13. — Normalized production factors for ^*Ti, ^"^Ni, and ^^Ni versus for adiabati c freeze outs from 



peak conditions CasA 1 — 3 in Table[SJ All points were drawn from a model for Cassiopia A (jMagkotsios et al 



20081 ) ■ Each central point represents a calculation that utilizes our recommended "semi-experimental" re- 
action rate for "^"Cafa, 7)' ^'*Ti and three choice s of ''^Ti(a,p)^^V reaction rate from our re-evaluation of 
the experimental effort of ISonzogni et alj (|2000l ): our "recommended" ^''Ti(a,p)'*^V rate (solid line, filled 
squares), and its upper and lower bound (denoted by filled triangles with dotted and dashed line styles, re- 
spectively). Error bars on each point reflect the minimum and maximum deviations of P44 due to variations 
of the six other '^°Ca(a, 7)"^^Ti reaction rates we consid ered. The dot-short-dash line shows '^''Ti synthesis 
using the original fits of Rauscher fc Thielemann ( 200Cll) for both '*°Ca(a, 7)^''Ti and ''^Ti(a,p)^''V. For the 
two NPF's P57 and P58 (shown as open and filled circles, respectively) an upper bound of 3 limits the allowed 
range of Ye- 
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Fig. 14. — Same as Figure fTSl but for pea k conditions 4 — 6 in Table [5] drawn from a rotating two-dimensional 
explosion model ( Magkotsios et al. 20081) . The production of ^''''^^Ni with respect to solar iron is nearly the 
same as in the model for CasA, but the production of ''^Ti is substantially lower in every case due to the 
lower entropy experienced in these expansions. 
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Fig. 15. — Evolution of select mass fractions versus temperature during freeze outs with identical peak 
temperatures (Tgp = 5.5) and electron fractions {Y^ = 0.4980), but different intial peak densities, for three 

points in our survey. The higher entropy expansions produce less ^^Ni and have a higher a-particle mass 
fraction at late times that facilitates the re-assembly of species below the iron group, including **°Ca and 
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Fig. 16. — Production factors (final decayed mass fractions divided by their abundances in the Sun) repre- 
senting the range of species produced in the a-rich freeze-out from peak conditions for point CasA-3 in Table 
[5] assuming our "recommended" reaction rates for ^°Ca(Q;, 7)'*^Ti and '*^Ti(a,p)^''V. This is for the specific 
case of l^e = 0.4980 and x = 1. This is a fairly typical nucleosynthesis pattern over the range of allowed 
for the a-rich freeze outs in the CasA and two-dimensional explosion models we consider. 
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Fig. 17. — Similar to Figurc[T6l this is the ratio ol the production factors for two scalings on the hydrodynamic 
time-scale, default (x = 1) and extended (x = 5). The species created in NSE are much less affected than 
the lighter isotopes which are reassembled for 2 < Tg < 4. 
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Fig. 18. — Temperature evolution in the CasA-3 expansion assuming two scalings on the hydrodynamic 
time-scale, default (x = 1) and extended (x = 5). Also shown is the ratio of the a-particle mass fraction. 
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Fig. 19. — Normalized production factors for Ti, ^'Ni, and ^^Ni versus Ye for adiabatic freeze outs from 
peak conditions CasA 1 — 3 in Table [5] but with an assumed x = 5 scaling of the hydrodynamic timescale. 
Results for the Ni isotopes are very similar to the default (x = 1) scaling case (Figure [T3|) . but due to the 
longer time spent at high temperature, the a-particle abundance at the onset of ^°Ca assembly was much 
lower 40%) than in the default scaling case, which facilitated a much lower final value for P44 in all cases 
studied. 



